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Abstract 


We  perform  finite-element  analysis  modeling  and  characterization  of  quasi-ballistic 
electrical  transport  in  semiconducting  carbon  nanotube  field-effect  transistors,  and  fit 
experimental  electrical  transport  data  from  both  suspended  and  on-substrate  single- 
walled  carbon  nanotube  transistors  fabricated  using  chemical  vapor  deposition.  Pre¬ 
vious  studies  have  focused  on  modeling  ballistic  transport  in  carbon  nanotube  field- 
effect  transistors,  but  have  ignored  the  spatial  dependence  of  the  resistance,  voltage, 
and  Fermi  energy.  These  spatial  variations  play  an  important  role  in  several  high- 
voltage  effects  that  are  particularly  important  in  the  quasi-ballistic  transport  regime 
where  most  current  or  near-term  devices  operate.  We  show  the  relationship  between 
device  geometry  and  pinch-off,  current  saturation,  and  channel  length  modulation  in 
the  quantum  capacitance  regime.  Output  resistance  is  found  to  increase  with  gate 
coupling  efficiency  with  a  power  law  behavior.  This  model  can  be  used  for  extraction 
of  device  properties  from  experimental  data  and  as  a  design  environment  tool. 
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1.  Introduction 


Single-walled  carbon  nanotube  (CNT)  field-effect  transistors  (FETs)  are  unique  devices  because  of 
their  intrinsically  one-dimensional  transport,  tunable  band  gap,  and  high  mobility.  They  also  offer 
potential  performance  advantages  such  as  high  linearity  in  analog  RF  amplifiers  [1];  high  surface-to- 
volume  ratio,  which  is  useful  for  chemical  sensors  [2];  high-frequency  mechanical  oscillations  for  RF 
NEMs  devices  [3];  and  high  cutoff  frequencies  in  the  THz  region  [4],  Several  research  groups  have 
used  theoretical  and  computational  methods  to  analyze  electrical  transport  in  CNT  FETs  [5],  [6],  most 
notably  Datta  et  al.  [7],  [8].  Several  software  tools  for  performing  these  simulations  are  available  on 
the  nanoFTUB.org  website,  such  as  the  FETToy  tool  [9],  Some  studies  have  focused  on  improving 
numerical  methods  [10],  [11],  whereas  other  researchers  have  focused  on  modeling  of  CNT  FETs  for 
process  design  guidance  [12]  and  circuit  design  [13].  Several  studies  have  emphasized  specific 
aspects  of  transport  in  CNT  FETs,  such  as  multiband  mobility  [14],  carrier  velocity  saturation  [15], 
Shottky  barrier  effects  [16],  and  electron-phonon  scattering  [17],  [18],  These  studies  and  tools,  how¬ 
ever,  all  assume  a  uniform  voltage,  resistance,  and  Fermi  energy  along  the  length  of  the  CNT.  For 
ballistic  transport,  this  assumption  is  not  valid  because  the  gate,  source,  and  drain  capacitances  vary 
spatially,  and  this  requires  a  spatial  solution  of  Poisson’s  equation  to  obtain  the  voltage  and  carrier 
number  density  distribution.  In  devices  that  experience  quasi-ballistic  and  diffusive  transport,  this 
approximation  precludes  modeling  of  pinch-off  and  channel  length  modulation,  effects  that  strongly 
influence  electrical  transport  in  FET  devices.  Here,  we  perform  numerical  modeling  of  electrical 
transport  in  CNT  FETs  within  the  Landauer  framework  while  explicitly  including  spatial  inhomoge¬ 
neity  by  calculating  the  voltage,  resistance,  and  Fermi  energy  at  discrete  points  along  the  CNT  chan¬ 
nel,  a  technique  known  as  finite-element  analysis.  Such  analysis  is  urgently  needed  in  order  to  exam¬ 
ine  the  dynamics  of  the  pinch-off  and  channel  length  modulation,  both  of  which  are  useful  to  gain  a 
complete  understanding  of  CNT  FET  operation  and  as  a  design  tool  for  device  processing  optimiza¬ 
tion  and  CNT-based  circuitry. 
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2.  Device  Fabrication  and  Experimental  Details 


In  order  to  validate  the  results  from  the  model  outlined  in  the  body  of  this  report,  experimental  data 
from  suspended  CNT  FETs  fabricated  using  chemical  vapor  deposition  (CVD)  are  compared  with 
output  from  the  model.  These  samples  are  fabricated  as  reported  previously  [19],  [20]  on  pre-defined 
platinum/tungsten  electrode  structures.  Figure  1  shows  scanning  electron  micrograph  (SEM)  images 
of  a  CNT  FET  device  fabricated  using  this  method.  The  CNT  is  grown  in  the  CVD  process  across  a 
trench  etched  in  the  underlying  substrate,  making  electrical  contact  on  both  sides  of  the  trench.  The 
trench  has  a  metallic  top  layer  that  serves  as  a  gate  electrode,  with  the  ambient  gas  or  vacuum  serving 
as  the  dielectric  between  the  CNT  and  the  gate.  The  yield  is  low  (~5%),  and  devices  with  CNTs  suc¬ 
cessfully  spanning  the  trench  are  selected  for  further  study.  The  CVD  process  is  the  last  step  in  fabri¬ 
cation,  ensuring  low  defect  concentration  in  the  resulting  CNT  channel  material. 
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Figure  1. 


CNT  device  geometry,  (a)  SEM  images  of  CNT  device  layout  and  (b) 
close-up  image  of  a  single  CNT  crossing  the  trench. 
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3.  Computational  Model 


In  the  Landauer  model,  electrical  current  flowing  in  a  ballistic  one-dimensional  conductor  is 

/  =  D(E)  ■  v g(E)[fR(E)-fL(E)]dE,  (1) 

where  vg(E)  is  the  group  velocity,  D(E )  is  the  density  of  states,  and  fR  and  fL  are  the  quasi-Fermi 
functions  for  the  right-  and  left-moving  electron  populations.  From  this  equation,  the  resistance  Rs  of 
a  diffusive  CNT  section  “s”  of  length  L  can  be  approximated  [14,  21]  as 

=  (2) 

where  R0  =  h/Ae2  is  the  quantum  resistance,  A  is  the  mean  free  path,  /  is  the  Fermi  function  for  the 
average  Fermi  energy  EF  =  (Ef  r  +  EF  L)/2  and  the  integral  is  taken  over  all  energies  except  the 

band  gap,  i.e.,  fbands  dE  =  dE  +  J^°  dE.  For  multiple  sections  of  one-dimensional  conductor  in 
series  with  one  another,  the  resistance  can  be  calculated  as 

R  =  *o  +  2/?s  +  2  Rc,  (3) 


where 


=  (4) 

is  :he  depleted  quantum  resistance  evaluated  at  the  point  of  minimum  doping  on  the  CNT,  and  Rc  is 
the  contact  resistance.  If  A  »  L  (ballistic  transport),  the  right-hand  side  of  (3)  reduces  to  the  depleted 
quantum  resistance  plus  the  contact  resistances.  The  mean  free  path  A  is,  is  given  by 

A  =  (Va  +  Va  )  (5) 

V  71 ac  ‘  /vop,emsJ 

where  Aac ,  the  mean  free  path  for  electron  scattering  by  acoustic  phonons,  is  given  by 

A„(E)=  (6) 

where  A^lfh  E  is  the  mean  free  path  at  high  energies  away  from  the  band  edge,  and  vF  is  the  carrier 
velocity  in  graphene.  Acoustic  scattering  is  approximated  as  elastic  as  in  [14].  A  more  complete 
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discussion  of  electrical  transport  modeling  with  inelastic  scattering  accounted  for  can  be  found  else¬ 
where  [22].  The  mean  free  path  is  related  to  the  scattering  rate  according  to  Aac(E]\  =  t(E)  * 

Vg(E).  Because  the  scattering  rate  increases  with  the  density  of  states  (Fermi’s  golden  rule)  near  the 
band  edge,  a  D(E)  a  vg~1{E').  A  full  derivation  of  this  effect  is  beyond  the  scope  of  this 

report,  but  the  divergence  of  the  scattering  rate  at  the  band  edge  where  the  group  velocity  goes  to  zero 
can  be  seen  elsewhere  in  the  carbon  nanotube  literature  [23-25].  The  combined  effect  gives  a  mean 
free  path  proportionate  to  vg.  The  energy  dependence  of  Aac  is  thermally  averaged,  but  still  varies  by 
a  factor  of  5  or  so  at  room  temperature  as  the  Fermi  energy  is  changed.  Electrons  in  one-dimensional 
systems  scatter  strongly  from  with  K-point  optical  phonons  [26].  The  mean  free  path  with  respect  to 
electron  scattering  through  optical  phonon  emission  Aopems  is  given  by 


A 


op,ems  ~ 


E  op  L 

~qV~ 


+  A 


op, min  ? 


(7) 


where  V  is  the  voltage  drop  across  the  segment,  Aop  min  is  the  scattering  length  for  electrons  with 
enough  energy  to  emit  an  optical  phonon,  and  Eop  =  0.1 6eV  is  the  optical  phonon  energy  of  the  zone¬ 
boundary  optical  phonons  [26].  The  voltage  drop  in  (7)  approximates  the  Pauli  exclusion  factor 
1  —  f{E  —  Eop).  This  relationship  was  discussed  by  Park  et  al.  [26],  and  can  be  justified  because 
kBT  «  Eop.  The  temperature  dependence  for  the  quantities  in  (5)-(7)  can  be  found  elsewhere  in  the 
literature  [27].  The  voltage  V(x )  and  EF(x )  at  given  section  along  the  CNT  can  be  calculared  by  solv¬ 
ing  Poisson’s  equation  [7]  to  obtain 

qEF(x)  =  ^ 


where  VG  is  the  gate  voltage;  =  CG  +  Cs  +  Cd  is  the  net  capacitance  from  the  section  of  the  CNT 
to  the  gate,  source,  and  drain  contacts;  and 

Q(Ef)  =  fjl  -  /(£)]  •  D(E)  dE  -  J“ /(£)  •  D(E)  dE  (9) 


is  the  linear  charge  density  on  the  CNT  section.  The  main  effect  of  the  source  and  drain  capacitances 
is  to  screen  the  influence  of  the  gate  near  the  contacts.  We  can  make  the  approximations  Cs(x )  = 

Csd0  —  and  Cd(x )  =  Csd0  ,  which  were  validated  by  an  electrostatic  finite-element  analysis  of  a 

X  xd~x 

standard  FET  device  geometry  in  Comsol  Multiphysics  software  package.  In  this  model,  Csd0  is  a 
characteristic  capacitance  per  unit  length  for  the  source-channel  and  drain-channel  capacitances,  xd  is 
the  source-drain  separation,  and  C0  is  the  gate  capacitance  per  unit  length.  Furthermore,  we  can  say 
that  Vs  =  Afi  and  Vd  =  A/j,  +  VT ,  where  A/x  is  the  surface  potential  of  the  contact  metal  relative  to  the 
CNT,  and  VT  is  the  voltage  drop  across  the  CNT. 

The  contact  capacitances  Cs  and  Cd  become  very  large  near  the  contact.  When  they  become  signifi¬ 
cantly  larger  than  both  the  quantum  capacitance  and  the  gate  capacitance,  (8)  dictates  that  EF  =  A/i, 
fixing  the  CNT  Fermi  energy  at  the  contact.  This  can  create  a  Schottky  barrier,  and  describes  the 
contact  doping  responsible  for  n/p- type  operation  in  CNT  FETs  [28].  This  energy  difference  is  not 
simply  the  work  function  difference  between  the  CNT  and  the  metal;  complex  surface  interactions 
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can  introduce  deviations  from  this  value  [29].  In  order  to  solve  this  equation,  EF (x)  and  Qix)  are 
solved  for  V (x)  numerically  with  (9)  by  use  of  a  lookup  table.  By  making  the  approximation  V ( x )  = 
—  VT  and  substituting  into  (8),  one  obtains 


qEF  — 


Q(gF)+AwAt(C5+Cd)+(VG-V(x))Cc 

CZ 


(10) 


This  approximation  assumes  a  linear  voltage  drop  along  the  CNT,  and  as  such,  fails  to  include  the 
effects  of  pinch-off  and  channel  length  modulation,  which  give  rise  to  finite  output  resistance  in  the 
saturation  region.  The  density  of  states  D(£)  is  given  by 


(ii) 


where  k(E)  is  the  wavenumber  in  the  hyperbolic  band  structure,  given  by 


k{E) 


(12) 


where  vF  —  840,000  m/s  is  the  Fermi  velocity  of  graphene,  and  Eg  is  the  band  gap. 

The  subthreshold  slope  S  is  related  to  the  gate  efficiency  according  to 

S  =  ——In  (10)  ,  (13) 

q  oc 


where  a  =  —  corresponds  to  the  gate  efficiency  of  the  segment  with  the  most  efficient  gate 

^Zmin 

capacitance  (the  least  contact  screening).  In  the  simplistic  model  for  the  contact  capacitance  outlined 
above,  Czmin  =  CG  +  4  Csd0.  The  voltage  drop  in  each  section  of  the  CNT  is  AT  =  IdRs,  where  Id  is 
the  drain  current.  In  order  to  perform  the  complete  calculation  for  the  electrical  conductivity  using  the 
equations  above,  the  device  parameters  needed  are  Eg,  Rc ,  X^lfhE ,  A/*,  C0,  oc ,  and  xd.  Given  experi¬ 
mental  data  from  a  CNT  FET  device,  it  is  possible  to  use  the  equations  above  in  a  fitting  exercise  to 
extract  these  parameters  from  the  device.  If  one  or  several  of  these  parameters  is  known  a  priori,  e.g., 
from  device  dimensions,  the  others  can  be  determined  with  greater  confidence.  Also,  since  the  rele¬ 
vant  physics  of  electrical  transport  are  being  considered,  these  equations  allow  for  realistic  modeling 
of  devices  with  known  parameters  in  a  design  environment. 

For  small  Vsd ,  the  current  and  spatial  voltage  profiles  of  the  CNT  are  calculated  for  a  given  Vsd  by 
iteratively  solving  the  above  equations.  This  approach  follows  the  general  scheme:  guess  V(x)  for  all 
x,  calculate  i?(x),  calculate  V(x)  again,  and  iterate  until  convergence  is  reached.  For  Vsd  high  enough 
to  induce  current  saturation,  however,  this  algorithm  is  problematic  because  for  each  iteration,  the 
calculated  drain  current  Id  changes,  which  results  in  a  change  in  V(x)  and  the  effective  gate  voltage, 
creating  an  unstable  feedback  loop.  Furthermore,  for  \Vsd  \  >  \  VG  —  Vth\,  parts  of  the  CNT  become 
completely  depleted  in  a  “super  pinch-off ?  condition,  and  the  simulation  produces  a  curve  that  exhib- 
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its  negative  differential  resistance  (NDR),  which  is  not  observed  in  field-effect  transistors  affected  by 
pinch-off  alone  (i.e.,  without  self  heating  [27]).  This  instability  is  mitigated  by  calculating  V  (x)  for  a 
fixed  Id  and  Vd,  and  using  an  error  minimization  algorithm  to  adjust  Id  until  V (xd)  =  Vd.  This  creates 
a  spatial  source-drain  causality,  where  the  voltage  of  any  given  element  is  only  affected  by  elements 
to  one  side.  This  procedure  is  more  computationally  intensive,  but  can  still  be  run  on  an  office  laptop 
computer. 
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4.  Results 


Figure  2  plots  experimental  data  (red  dots)  taken  from  an  isolated,  suspended,  semiconducting  CNT 
FET  showing  drain  current  plotted  as  a  function  of  gate  voltage.  The  device  is  modeled  as  a  p- type, 
with  high  on-current  at  negative  gate  voltages,  an  on-off  ratio  of  approximately  1  x  106,  and  a  sub¬ 
threshold  slope  of~100  mV/decade.  Plotted  along  with  the  experimental  data  is  a  curve  calculated 
using  Landauer  transport  theory  (blue  line).  In  order  to  fit  the  data,  the  simulation  parameters  are: 
gate  capacitance  CG  =  lOpF/m,  gate  coupling  efficiency  a  =  62%,  high-energy  mean  free  path  A^hE 
=  1100  nm,  band  gap  Eg  =  0.85  eV,  and  contact  resistance  Rc  =  60  kft.  The  CNT  channel  length  was 
set  at  xd  =  1  pm,  in  accordance  with  the  actual  device  geometry.  For  VG  >  IT,  the  theoretical  model 
does  not  agree  well  with  the  experimental  data.  Plere,  the  CNT  FET  undergoes  bipolar  transport 
(diode  behavior,  with  p-typz  contacts  and  n- type  channel),  which  is  not  considered  by  the  model.  A 
further  discussion  of  electrical  transport  can  be  found  elsewhere  in  the  literature  [30,  31]. 

Figure  3  shows  a  curve  manifold  containing  experimental  Id  versus  Vsd  data  measured  from  a  sub¬ 
strate-bound  CNT  FET  from  Zhang  et  al  [32]  (reprinted  with  permission),  fabricated  with  Hf02  gate 
dielectric  and  Pd  contact  electrodes.  Zhang  et  al  state  that  the  channel  length  was  4  pm,  the  diameter 
2  nm,  and  the  gate  capacitance  150  pF/m.  The  I-V  curves  were  measured  at  a  variety  of  gate  voltages 
from  -2.6V  to  -0.6V,  and  are  plotted  along  with  simulation  results  from  the  theoretical  model  out¬ 
lined  here.  The  device  is  p- type,  with  larger  conductance  at  negative  gate  voltages.  The  behavior  is 
ohmic  at  low  voltages  before  transitioning  into  the  saturation  region.  The  simulation  parameters  were 


Ks(V) 

Figure  2.  Low- voltage  electrical  transport  data.  Drain  current  versus  gate  voltage 
for  a  semiconducting  CNT.  Experimental  data  (red  circles)  and  theoreti¬ 
cal  transport  model  (blue  line)  agree  nicely  for  VG  <  1 V. 
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Figure  3.  High-voltage  electrical  transport  data.  Experimental  (open  cir¬ 
cles)  and  theoretical  (blue  lines)  drain  current  plotted  versus 
source-drain  voltage  for  various  gate  voltages.  The  experimental 
data  is  from  Zhang  et  al  [32]  (reprinted  with  permission). 

manually  adjusted  until  a  reasonable  fit  was  achieved.  The  simulation  parameters  are:  gate  capaci¬ 
tance  CG  ~  120  pF/m,  gate  coupling  efficiency  a  -  99%,  high-energy  mean  free  path  =  580 

nm,  band  gap  Eg  =  0.42  eV,  contact  resistance  Rc  =  3.96  k£2,  values  which  agree  reasonably  well 
with  the  experimental  device  parameters.  The  CNT  channel  length  was  set  at  xd  =  4  jam,  in  accord¬ 
ance  with  the  value  reported  by  Zhang  et  al.  The  fitting  parameters  introduced  here  are  the  minimum 
number  of  parameters  necessary  to  capture  the  physics  of  channel  length  modulation  in  a  CNT  FET. 
Similar  current  saturation  plots  have  been  reported  in  the  literature  [32-37].  Due  to  the  difficulty  of 
fitting  so  many  different  data  curves  through  trial  and  error,  there  is  some  deviation  between  the  theo¬ 
retical  model  and  the  experimental  data.  Most  of  the  deviation  is  concentrated  at  large  negative  Vq 
and  large  VS(/,  when  the  device  is  deep  in  the  “ON”  state.  Despite  these  small  deviations,  the  theoreti¬ 
cal  model  fits  the  experimental  data  reasonably  well.  Computer-automated  error  minimization  algo¬ 
rithms  could  reduce  this  error;  however,  such  investigation  was  beyond  the  scope  of  this  work. 

Figure  4  shows  the  linear  resistivity  along  the  length  of  the  CNT  at  selected  Vsd ,  taken  at  a  gate  volt¬ 
age  ofVG  =  —1  V,  for  the  same  substrate-bound  CNT  FET  of  Figure  3.  The  curves  illustrate  the  pinch- 
off  effect,  which  is  responsible  for  current  saturation  at  high  source-drain  bias  voltages.  Near  the 
drain  electrode,  the  voltage  drop  along  the  CNT  causes  the  effective  gate  voltage  to  be  lowered, 
reducing  the  charge  density  in  the  pinch-off  region,  which  becomes  progressively  more  and  more 
depleted  of  carriers.  As  a  result,  the  resistivity  in  this  region  is  large.  The  charge  carrier  densities  in 
sections  of  the  CNT  channel  to  the  left  of  the  pinch-off  region,  however,  remain  relatively  unchanged. 

Figure  5  shows  drain  current  saturation  curves  for  the  same  CNT  FET  with  different  values  for  gate 
efficiency  a.  For  high  gate  efficiency  (low  channel-source  and  channel-drain  capacitances),  the  out¬ 
put  resistance  is  high  with  a  finite  output  resistance  in  the  current  saturation  region  due  to  channel 
length  modulation  caused  by  channel-drain  capacitive  coupling.  For  low  gate  efficiencies,  however, 
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Figure  4.  Pinch-off  in  a  CNT  FET.  Resistivity  versus  spatial  position  along  the 
length  of  a  CNT  FET  during  development  of  the  pinch-off  condition. 
Inset  shows  drain  current  versus  source-drain  voltage  for  this  calculation. 


^  (V) 

Figure  5.  Channel  length  modulation  in  a  CNT  FET.  Drain  current  versus 
source-drain  voltage  at  gate  efficiencies  a  of  99%,  98%,  and 
95%,  showing  decreased  output  resistance  for  lower  gate  cou¬ 
pling  efficiencies  caused  by  channel  length  modulation. 

the  output  resistance  is  low.  From  the  equations  above,  the  gate  efficiency  is  described  by 


ii 


The  IV  curves  shown  in  Figure  5  were  calculated  for  many  different  values  of  CG ,  a,  and  Vq\  and  the 
slope  of  the  current  saturation  region  in  the  resulting  IV  curves  was  fit  to  determine  the  output 
resistance  as  a  function  of  the  different  variables  in  order  to  examine  the  relationship  between  device 
geometry  and  pinch-off,  current  saturation,  and  channel  length  modulation.  The  range  of  gate  capaci¬ 
tance  modeled  in  the  simulation  includes  values  from  15  pF/m  to  480  pF/m.  This  range  was  chosen 
deliberately  to  include  values  both  above  and  below  the  theoretical  quantum  capacitance  for  a  CNT  of 
96  pF/m  [38,  39].  Figure  6  shows  Routput  plotted  versus  1  -  a  for  different  gate  voltages,  with  CG  = 

15  pF/m  and  gate  voltage  range  -2.6V  <  VG  <  -0.6V.  The  threshold  voltage  Vth  ~  -0.6  V.  R0utput 
displays  a  power  law  dependence  on  1  -  a  according  to  Routput  (1  “  a)p  with  P  =  “0.93  at  VG  ~ 
Vth  and  p  =  -0.69  at  VG  =  -2.6V.  For  higher  gate  capacitances  the  value  of  p  varies  more  widely.  At 
CG  =  480  pF/m,  p  =  -1.0  at  VG  ~  Vth  and  p  =  -0.38  at  VG  =  -2.6V. 

From  the  numerical  simulations,  it  is  also  possible  to  calculate  the  transconductance  gm  for  different 
transistor  configurations  in  the  current  saturation  regime.  Figure  7  shows  the  calculated  gm  versus  VG 
for  the  same  CNT  FET  with  several  different  values  of  Cg,  a  gate  efficiency  a  =  96%,  and  Vsa  =  -1.5 
V.  The  transconductance  peaks  at  different  voltages  for  the  different  gate  capacitances,  with  trans¬ 
conductance  actually  decreasing  with  increasing  gate  capacitance  for  VG  <  -1.6V.  Possible  causes  for 
this  could  be  saturation  of  the  CNT  channel  for  high  carrier  concentrations,  or  optical  phonon  emission 
onset  at  the  high  drain  currents  that  accompany  high  carrier  concentrations.  gm  was  also  characterized 
with  respect  to  the  gate  efficiency  a;  however,  it  was  found  to  vary  very  little  with  changes  in  a  . 


0.01  0.1 
1-a  (unitless) 


Figure  6.  Output  resistance  as  a  function  of  1  -  a  for  different  gate  volt¬ 
ages  (Vg-  - 2.6  V  to  -0.6  V).  Routput  has  a  power  law  depend¬ 
ence  on  1  -  a  ,  with  Routput  increasing  as  a  ->  1 . 
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Figure  7.  Transconductance  gm  versus  gate  voltage  for  several  different 
gate  capacitance  values  at  a  Vsd=-\.5  V  and  a  =  96%. 


13 


5.  Conclusion 


In  conclusion,  we  have  used  the  Landauer  model  used  to  calculate  electrical  transport  in  semicon¬ 
ducting  carbon  nanotube  field-effect  transistors  at  high  source-drain  voltages  and  for  arbitrary  gate 
voltages  in  the  quantum  capacitance  limit.  The  model  is  used  to  fit  experimental  electrical  transport 
data  from  suspended  and  substrate-bound  single-nanotube  CNT  FETs.  The  model  successfully  repro¬ 
duces  the  subthreshold  slope,  on-off  ratio,  superthreshold  roll-off,  and  current  saturation  at  high  volt¬ 
age.  The  model  captures  the  physics  of  the  channel  length  modulation  effect,  which  results  in  finite 
output  resistance  in  FETs.  This  model  is  a  useful  tool  in  characterizing,  predicting,  and  optimizing 
electrical  performance  of  carbon  nanotube  field-effect  transistors  for  both  digital  and  analog 
applications. 
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